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Abstract 

We present a rigorous next-to-leading order analysis of J /if/ yield and polarization within the factorization theorem 
of nonrelativistic QCD (NRQCD). To the orders considered, this framework depends on three free parameters, the 
color-octet long-distance matrix elements. We extract their values in a global fit to inclusive J/ ifi production data from 
various hadroproduction, photoproduction, two-photon scattering and electron-positron annihilation experiments. We 
show that this fit is constrained and stable and describes all data sufficiently well. We then make predictions for J '/if/ 
polarization in photo- and hadroproduction and compare them to the currently available data. As for photoproduction, 
HERA data is not precise enough to draw definite conclusions. But as for hadroproduction, CDF data measured at 
Tevatron run II is in strong conflict with NRQCD predictions. With early ALICE data being however compatible 
with NRQCD, the future, more precise polarization measurements at the LHC will thus have the potential to clearly 
confirm or dismiss LDME universality. 

Keywords: Lepton-Nucleon Scattering, Nucleon-Nucleon Scattering, Lepton-Lepton Interactions, Heavy quarkonia, 
Nonrelativistic QCD 



1. Introduction 

Heavy quarkonia are bound states of a heavy quark 
and its antiquark. There are charmonia and bottomonia. 
The charmonium J/if/ is the most extensively studied 
quarkonium state, because of its experimentally clean 
signature due to the large branching fraction of its lep- 
tonic decay modes. Almost 40 years since its discovery, 
the underlying mechanisms governing heavy quarko- 
nium production and decay are however still not under- 
stood and subject to dispute. According to the factoriza- 
tion theorem of nonrelativistic QCD (NRQCD) fl], the 
cross section to produce a heavy quarkonium H factor- 
izes according to 

o-(ab -> H + X) = ^ cr(ab -> cc[n] + X)(0 H [n]), 

n 

where the cr(ab — > cc[n] + X) are perturbatively cal- 
culated short distance cross sections describing the pro- 
duction of a heavy quark pair (here cc) in an intermedi- 
ate Fock state n, which does not have to be color neu- 



tral. The (0 H [n]) are nonperturbative long distance ma- 
trix elements (LDMEs), which currently have to be ex- 
tracted from experiment and which describe the tran- 
sition of that intermediate cc state into the physical H 
via soft gluon radiation. NRQCD predicts each of the 
LDMEs to scale with a definite power of the relative 
heavy quark velocity v, which serves as an additional 
expansion parameter besides a s : In case of H = J/ if/, 
the leading order contribution in the v expansion stems 
from n - S l and equals the traditional color sin- 
glet model (CSM) prediction, while the leading rela- 
tivistic corrections are made up by the l S^\ 3 S' SI , and 
3 Py ] states. The upper index "8" stands for color octet 
(CO), and these contributions are usually just called the 
color octet contributions. The aim of our work is to test 
the universality of the LDMEs and thereby to challenge 
NRQCD factorization. 
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2. History of NLO J /if/ production calculations 

Next-to-leading order (NLO) corrections to inclu- 
sive J /iff production cross sections have been calcu- 
lated within the CSM for direct photoproduction (2), for 
hadroproduction [3| and for electron-positron annihila- 
tion El. 

The calculation of NLO corrections to the short dis- 
tance cross sections of the intermediate CO states, espe- 
cially to the 3 f^ 8] states, have on the other hand proven 
to be very challenging. But up to now they have been 
calculated for all relevant collision processes as well. 
The 2 — > 1 processes for photo- and hadroproduction 
have already been calculated in 1998 [5 1. Inclusive pro- 
duction in direct two-photon collisions followed in 2005 
(6|, in direct photoproduction [7] and electron-positron 
scattering neglecting the small 3 5j 81 contribution [8l in 
2009. The hadroproduction calculations were still 
missing the ] contributions. Full calculations in- 
volving all CO states followed in 2010 with two in- 
dependent works IflOl ITT1 . The missing pieces of sin- 
gle and double resolved two photon scattering, resolved 
photoproduction and the 3 S j contributions of electron- 
positron scattering were finally presented in 2011 fl2l . 

Polarized NLO J /if/ production cross sections have 
been evaluated within the CSM, for direct photoproduc- 
tion |[T3l and hadroproduction [14|, and for the and 
3 5 [ 1 S] intermediate states in hadroproduction |9|. Re- 
cently, polarized NLO cross sections including all CO 
contributions have been calculated, namely for direct 
photoproduction at HERA 1 15 1 and for hadroproduction 
at Tevatron and the LHC Q6). 

3. J /if/ yield: A global fit of the CO LDMEs 

In IPT21 we have published a global NLO fit of the 
three CO LDMEs to 194 data points of inclusive unpo- 
larized J /if/ production from 10 different experiments, 
see figure [TJ and table [JJ for the fit results. These ex- 
periments include data from photoproduction at HERA, 
hadroproduction at RHIC, Tevatron and the LHC and 
additionally data from two-photon collisions measured 
at LEP and electron-positron collisions at KEKB. Most 
of the data fitted to is prompt, while none of our short 
distance cross sections does include feed-down contri- 
butions from higher charmonium states. Therefore, in 
table [JJ we list two sets of CO LDMEs: In set A, we 
ignore these contributions, this is our default fit, while 
in set B we take care of them by subtracting estimated 
feed-down contributions from the con-esponding data 
prior to fitting. The values of set A are used for all 
curves of figures [JJ [2] [3] and [4] In figure [2] and [3] some 



diagrams of figure [JJ are decomposed into contributions 
of the different intermediate Fock states as well as into 
direct and resolved photon contributions. The global 
fit shows that at NLO all considered processes except 
perhaps the two-photon collisions can be described well 
with a single set of CO LDMEs. Furthermore, the CO 
LDMEs thus obtained do indeed exhibit the scaling be- 
havior predicted by the scaling rules of NRQCD HI: 
They are roughly of 0(v 4 ) relative to {O'^^S [ 1] )>, with 
v 2 w 0.2. The color-singlet (CS) contributions alone 
are on the other hand shown to fall clearly short of the 
data everywhere except for the BELLE total <? + e~ cross 
section. 

At this point let us look at the different production 
processes in more detail. Figures [T^-o show J/if/ trans- 
verse momentum pj distributions of various hadropro- 
duction measurements. Nearly all data points lie within 
the NRQCD NLO band, which estimates corrections 
due to even higher orders in a s . From a comparison 
of the LO and NLO curves the a s expansion seems to 
converge much more rapidly when the CO contributions 
are included than in the CSM, where the slope drasti- 
cally changes when going from LO to NLO. At higher 
values of transverse momentum than considered in the 
fit, resummations of large logarithms \og(p T /Mj^) ap- 
pearing in the short distance cross section will become 
necessary, as can be seen from our comparison to re- 
cent ATLAS data in figure [4] Such resummations are 
a standard tool to enhance the applicability of fixed- 
order calculations in multi-scale problems, see figure 8 
of [30 1 as an example. In the case of quarkonium pro- 
duction, this might best be achieved within the frame- 
work outlined in GTJ. In figure |2ji, the NLO CS+CO 
cross section of plot[JJ;, as an example, is decomposed 
into the various contributions of the different interme- 
diate cc Fock states. These are the short distance cross 
sections already multiplied by their respective LDMEs. 
The fact that the overall 3 fJ / 8] contributions are negative 
for pj J$ 7 GeV (here the short distance cross section 
is positive, but the LDME is negative, for pj ^ 7 GeV 
both are negative) is not worrying since individual con- 
tributions are unphysical, only the overall cross section 
has to be positive. 

The global fit does also describe the photoproduction 
at HERA well. In the photoproduction limit, the incom- 
ing electron or positron interacts with the proton via a 
quasi-real bremsstrahlung photon. Figures [TJ— z show 
distributions in the J /if/ transverse momentum pj, in the 
photon-proton invariant mass W and in the inelasticity 
variable z, which in the proton rest frame is the fraction 
of the photon energy taken over by the J/if/. Here, even 
the z distribution is now much better described than in 
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Figure 1: Plots a-t: Results of the global fit (T3) compared to ALICE (T7), ATLAS (TS), Belle (T5), CDF (30]|IT], CMS (22), DELPHI (23), LHCb 
1 24 1 , PHENIX [25], and ZEUS |26] data. The blue bands are the CSM predictions, the yellow bands include the CO contributions. The bands are 
constructed by variation of the renormalization, factorization and NRQCD scales. 
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Figure 1: Plots u-z (continuation): Results of the global fit 1 12 ] compared to HI [27 28] data. The blue bands are the CSM predictions, the yellow 
bands include the CO contributions. The bands are constructed by variation of the renormalization, factorization and NRQCD scales. 




Figure 2: Decomposition of the NLO CS+CO curves from figures[T];, q, r, and t into the contributions of individual intermediate cc states. The line 
coding is the same for all four graphs. Red curves mean negative values. Please note that these curves are the short distance cross sections already 
multiplied by the corresponding LDMEs. 




Figure 3: Figure a shows a decomposition of the NLO CS+CO curve of the two-photon-scattering results from plot^j into direct, single- and 
double-resolved photon contributions. Figures b and c show decompositions of the NLO CS+CO curves of the pj and z distributions in photopro- 
duction shown in figuresfTj' and t into direct and resolved contributions. 
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Set A: Do not mind feed-downs 


Set B: Subtract feed-downs first 


(0 J/ >"( l S l * s )) 

(0«( 3 r')) 
(0 Jl f( 3 P [ o ] )) 


(4.97 + 0.44) x 10" 2 GeV 3 
(2.24 ± 0.59) x 10~ 3 GeV 3 
(-1.61 ±0.20)x 10~ 2 GeV 5 


(3.04 ± 0.35) x 10~ 2 GeV 3 
(1.68 ± 0.46) x 10~ 3 GeV 3 
(-9.08+ 1.61) x 10" 3 GeV 5 



Table 1: Results of global fit ['12'] for the J/ip CO LDMEs. Set A corresponds to the main fit shown in figure [T] In set B, estimated feed-down 
contributions from higher charmonium states were subtracted from the prompt data prior to fitting (hadroproduction: 36%, photoproduction: 15%, 
yy scattering: 9%, e + e~ annihilation: 26%). 
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Figure 4: Predictions for the ATLAS measurement 1291 . Note that 
these data are not part of our global fit, since they became public after 
our global fit was finished. The bands are again constructed by varia- 
tion of the renormalization, factorization and NRQCD scales. At very 
high pr, it will be necessary to resum large logarithms log(pj /npj^). 
For instance, at p T = 40 GeV, a s logG^/Af^J a 0.7. 



previous analyses, for two reasons: First, this analysis 
is the first one to include resolved photon contributions, 
in which the photon in turn interacts with the proton via 
its hadronic content. Resolved photoproduction domi- 
nates the cross section below z ~ 0.25, see figure [3J;. 
And secondly, the cross section is at high z now much 
better described than in the older Born analyses, which 
predicted a steep rise in the cross section not found in 
the data. The reason is a strong cancellation between 
the ! 5q 8] and 3 f', 8] contributions, due to the negative 
value of (O^^CP 1 ^)}, as shown in fi eureEti. We further 
note that the HI HERA 2 data (28] (fig.gjf-z) lie sys- 
tematically below the data measured by HI at HERA 1 
[27] (fig.[2ji-w) and by the ZEUS collaboration both at 
HERA 1 [26) (fig.g-t) and HERA 2 l32l . 

Figure [T]j shows the transverse momentum distribu- 
tion for two-photon scattering as measured by the DEL- 
PHI collaboration at LEP Il23l . In this global analysis 
lfT2ll . for the first time also NLO contributions due to sin- 
gle and double resolved photons have been calculated. 
The single resolved contributions are in fact contribut- 
ing up to 99% to the overall cross section, see figure [3^. 



The fact that the DELPHI data overshoots the NRQCD 
prediction is not worrying since the experimental errors 
are huge with just 16 events entering the data of fig- 
ure [T}j. 

As for //t// production in electron-positron annihila- 
tion, the production in association with another charmo- 
nium or open cc state is dominating the inclusive cross 
section lfl9l . Since these contributions are not included 
in our calculation, in figure [TJ) we compare our results 
for the total production cross section at KEKB with 
the measurement lfT9ll performed by Belle, in which the 
double charmonium and J /iff + cc contributions are ex- 
plicitly subtracted from the inclusive cross section. We 
see here that both the NLO CS and the NLO CS+CO re- 
sults are compatible with the data. However, two things 
have to be noted here: First, in all our calculations, the 
partonic final states of the LO amplitudes consist of //<// 
plus one light quark. But within the CSM, the Born pro- 
cess is e + e~ — > J/t//+gg, which in our terminology is an 
NLO processes. We do not consider NLO corrections to 
e + e~ — > J/ift + gg for reasons of consistency, although 
they have been shown to further increase the CSM pre- 
diction |4|. The second caveat is that in the Belle anal- 
ysis lfT9l only events with more than 4 charged tracks 
have been measured. The correction factor making up 
for the missing events could not be determined, there- 
fore the actual cross section is larger than the presented 
one by an unknown factor. 

3.1. Dependence on low px cuts 

In our default global fit, we impose a low-p T cut on 
all hadroproduction data of 3 GeV, since data with lower 
Pt can not be described in a fixed-order treatment due 
to soft gluon radiation: The data exhibits a flattening, 
which can not be successfully described. For a simi- 
lar reason we consider only two-photon scattering data 
with p T > 1 GeV. In this section we show the stability of 
our global fit with respect to varying the low-pj- cuts on 
the data. In that way we also show that our global fit is 
constrained, while a fit to hadroproduction data alone is 
not. We also successfully compare to the Tevatron-only 
fit published in |[T0ll33l 
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p T > 1 GeV 


p T > 2 GeV 


p T > 3 GeV 
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148 points 


134 points 


119 points 


86 points 


60 points 


(O^Csp) [io- 
(o^Csf 1 )) [io- 


' 2 GeV 3 ] 


5.68 ± 0.37 


4.25 ± 0.43 
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3.91 ±0.51 


3 GeV 3 ] 


0.90 ± 0.50 


2.94 ± 0.58 
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2.23 ± 0.62 
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2 GeV 5 ] 


-2.23 + 0.17 


-1.38 ±0.20 


-1.61 ±0.20 


-1.59 ±0.22 


-1.16 ±0.23 



Table 2: Global fit results for different low-pj- cuts on the hadroproduction data. Besides the number of hadroproduction data points listed, the fits 
include 74 data points of photoproduction and two-photon scattering, and 1 data point of e + e~~ annihilation. The fit with pj > 3 GeV is our default 
fit. One can see that the fit results are practically independent of the low-p^ cut, unless the cut is taken below 2 GeV. 



p T > 1 GeV p T > 2 GeV p T > 3 GeV p T > 5 GeV p T > 7 GeV 

Photoproduction and 77 data left 74 points 30 points 15 points 5 points 1 points 

{O^CS 1 ^)) [10"2 GeV 3 ] 4.97 ± 0.44 5.10 ± 0.92 4.05 ± 1.17 5.44 ± 1.27 9.56 ± 1.59 

(tf^fSp) [10~ 3 GeV 3 ] 2.24 ±0.59 2.11 ± 1.22 3.52 ± 1.56 1.73 ± 1.68 -3.66 ±2.09 

<0W(3p[8])} [10 -2 GeV 5] -1.61 ±0.20 -1.58 ±0.48 -0.97 ±0.63 -1.63 ±0.68 -3.73 ± 0.83 

Table 3: Global fit results for different low-pr cuts on the photoproduction and two-photon scattering data. Besides the number of photoproduction 
and two-photon scattering data points listed, the fits include 1 19 data points of hadroproduction, and 1 data point of e + e~ annihilation. The fit with 
Pt > 1 GeV is our default fit. The necessity of a critical amount of photoproduction data to stabilize the fit is clearly seen. In the region of a low- pr 
cut between 1 GeV and 5 GeV, where there is enough photoproduction data included, the fit results are practically independent of that cut. 
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148 points 



p T > 2 GeV 
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p T > 3 GeV 
119 points 



p T > 5 GeV 
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p T > 7 GeV 
60 points 



(Cy^CSf)) [10~ 2 GeV 3 ] 
((y/^S™)) [10~ 3 GeV 3 ] 
(0W(3pjM)) [iQ-2 GeV 5 ] 



8.54 ±0.52 16.85 ± 1.23 11.02 ± 1.67 1.68 ±2.20 2.18 ±2.56 
-2.66 ±0.69 -13.36 ± 1.60 -5.56 ±2.19 8.75 ± 2.98 10.34 ± 3.55 
-3.63 ±0.23 -7.70 ±0.61 -4.46 ±0.87 2.20 ± 1.23 3.50 ±1.50 



M [10~ 2 GeV J ] 
Mi [10~ 3 GeV 3 ] 



2.25: 
6.37: 



:0.12 
:0.19 



3.51 : 
5.80: 



0.19 
0.19 



3.29 ± 0.20 
5.54 ± 0.20 



5.50 ± 0.29 
3.27 ± 0.29 



8.24 ± 0.58 
1.63 ±0.43 



Table 4: Fits restricted to only hadroproduction data: Fit results for different \ow-pr cuts. The fits include only the listed number of hadroproduc- 
tion data points. The fits are underconstrained, and the fit results depend strongly on the low-pj- cut chosen. For comparison, also the values for 
the linear combinations Mq and M\ used in references 1 10 33 ] are given. 
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In table [2] we list our global fit results for different 
values of the low-pr hadroproduction cut. We see that 
the fit results are almost independent of this cut. The 
values for the fitted LDMEs vary little for low-p T cuts 
down to about 2 GeV. In table [3] we then show results 
of our global fit with different cuts on the photoproduc- 
tion and two-photon scattering data. We see that the fit 
is stable against varying this cut in the region between 
1 GeV and 5 GeV. However, with even higher low-p T 
cuts, too little photoproduction data is left, so the fit be- 
comes unstable. 

The necessity to include data from different produc- 
tion mechanisms becomes even more obvious when 
we vary the low-pr cut in fits to hadroproduction data 
alone. In table [4] we show the corresponding results 
when we restrict our global fit to hadroproduction data. 
The fit results show a very strong dependence on the 
low- pr cut, and the 3 -parameter fit to only hadropro- 
duction data is clearly underconstrained. Therefore, in 
|fT0ll33l , where a fit to CDF Tevatron run II data alone 
was performed, the authors fit the two linear combi- 
nations M = (O^CS 1 ^)) + 3.9(0 / ^( 3 J P[ ) 8] )>/w2 and 
M 1 = {0 JI *CS { * ] )) - 056{O JI *CP { ^))lm 2 c . For com- 
parison we list our results for Mo and M\ in table [4] as 
well. We note that the results of our fit to only hadropro- 
duction data with pj > 7 GeV agree well with their 
corresponding results 

M = (8.54 ± 1.02) x 1(T 2 GeV 3 , and 
Mi = (1.67 ± 1.05) x 10~ 3 GeV 3 , 

which are listed in equation (18) of |33|. 

4. Predictions for polarization observables 

In order to scrutinize the universality of the LDMEs 
it is however not enough to consider only the inclu- 
sive J/if/ yield. Therefore we now apply the set of CO 
LDMEs extracted in the last section to make predic- 
tions for J /iff polarization observables. More specifi- 
cally, we use set "B" of table [T] The results for photo- 
and hadroproduction are summarized in this section. 

Measuring the polarization of J /if/ means measuring 
the angular distribution of the two leptons, by which the 
J /iff is tagged. This distribution is parameterized via 

W(G, <p) oc l+A cos 2 O + Aj sin 2 6 cos(20) 
+ Ag,p sin(26>) cos <f>, 

where 6 and <p are the polar and azimuthal angles of the 
p + or e + in the J /if/ rest frame. This definition does 
of course depend on the choice of the coordinate sys- 
tem axes. Among the frequently used coordinate frames 



are the helicity frame, in which the polar axis points in 
the direction of -(p y + p p ), the Collins-Soper frame, 
in which it points to p y /\p y \ - p P /\p P \ and the target 
frame, in which it points to -p p . Here, p p is the three 
momentum of a colliding proton, and p y the three mo- 
mentum of the colliding photon in case of photoproduc- 
tion, and of the second proton or the antiproton in case 
of hadroproduction. Ag = corresponds to unpolar- 
ized J /if/, whereas A g = +1 (-1) stands for fully trans- 
versely (longitudinally) polarized J /if/. Unfortunately, 
there are different naming conventions for the polariza- 
tion parameters. In the photoproduction literature, the 
parameters A = Ag, fi = Ao^ and v = 2A^ are widely 
used, and on top of that Ag is often called a. 

On the theoretical side, we calculate the parameters 
Ag, A,/, and A gt/ , via 

_ dan -da-QQ _ &yi 

do~\ i + do~oo ' ^ d<T\ i + da-QQ ' 
^ V2Reafcr 10 

^60 _ 3 j ' 

acT\\ + clcr go 

where dcr-,j are the differential J/ if/ production cross sec- 
tions, calculated using the NRQCD factorization, but 
keeping the spin of the intermediate cc[n] pair fixed in- 
stead of summing over it. The spin polarization vectors 
e*(i) in the amplitude and e(j) in the complex conju- 
gated amplitude are instead replaced by their explicit 
expressions, as for example derived in [34|. In case of 
the Fock state n = 3 P^f \ for which S — L — 1, the or- 
bital angular momentum is, however, still summed over. 
For the spin-zero Fock state n = l S Q & \ dcroo and dtr\\ 
are each set to one third of the unpolarized cross section 
and d(T\-\ — do~ w — 0. 

4.1. Results for photoproduction 

Our results for direct photoproduction [ 15 1 are shown 
in figure [5] We compare our predictions for the pa- 
rameters A and v as functions of p T and the inelasticity 
variable z with data measured by the HI [28 1 collabora- 
tion in the helicity and Collins-Soper frames and by the 
ZEUS collaboration [35] in the target frame. 

We note that the direct photoproduction mechanism 
is the dominant mechanism only in the intermediate z 
region 0.3 z J$ 0.9. At lower z, resolved photopro- 
duction starts to dominate (see figure [3};), and at higher 
z diffractive J /if/ production overwhelms the cross sec- 
tion. Furthermore, close to the kinematical endpoint 
region z ~ 1, the NRQCD v expansion is expected to 
break down. Although the ZEUS pj distribution mea- 
surement does not impose an upper z cut, we therefore 
nevertheless integrate only up to z = 0.95 in figure [5};. 
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Figure 5: Polarization parameters A and v for direct photoproduction at HERA using CO LDME set B of table[T] compared to HI (28] and ZEUS 
1 35 1 data. Blue bands: Uncertainties of NLO CS curve due to scale variations. Yellow bands: Uncertainties of NLO CS+CO curve due to scale 
variations and uncertainties of the CO LDMEs. From 1151 . 
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We note that the strong tendency towards transversely 
polarized J /if/ in plot [5]; could possibly be related to 
exactly that feature, since diffractively produced vector 
mesons are indeed predicted to be strongly transversely 
polarized in the endpoint region ll36ll . 

Unfortunately, the HI (28] and ZEUS £35) data do 
not yet allow to distinguish the production mechanisms 
clearly, although the overall^ 2 value of all the HI and 
ZEUS data in figure [5] with regard to the default NLO 
predictions is reduced by more than 50% as the CO con- 
tributions are included. But kinematical regions can be 
identified, in which a clear distinction could be possible 
in more precise experiments at a future ep collider, like 
the LHeC: At higher p T , NRQCD predicts the J/if/ to 
be largely unpolarized in contrast to the CSM predic- 
tions. In the z distributions, however, the scale uncer- 
tainties are sizeable and the error bands of the CSM and 
NRQCD predictions largely overlap. 

4.2. Results for hadroproduction 

Our results for hadroproduction IT6l are shown in fig- 
ure [6] We compare our predictions for the parameters 
Ag and A^ as functions of pj in the helicity and Collins- 
Soper frames with the measurements by CDF 11371 [381 
and ALICE (39). Analogously to the unpolarized case 
we do not consider the p T range below 3 GeV, where 
nonperturbative effects are likely to dominate. In fig- 
ure|7]we decompose the contributions to dcroo, dan and 
dcr\ _i used for the predictions of the ALICE measure- 
ments into the contributions of the various intermedi- 
ate Fock states. The unpolarized cross section is recov- 
ered as dcroo + 2d<j\ \ . Our predictions do not yet in- 
clude feed-down contributions from higher charmonium 
states, whereas the CDF data are prompt, and the AL- 
ICE data even non-prompt, but the leading-order (LO) 
NRQCD analysis 11401 has shown the impact of Xcj and 
if/' feed-downs to be rather small. In the helicity frame, 
the CSM predicts strongly longitudinally polarized J /if/ 
at NLO, while NRQCD predicts a strong transverse po- 
larization. In the Collins-Soper frame the situation is 
inverted. The CDF measurement at Tevatron Run II 
l3lfl . which is partially in disagreement with the mea- 
surement at run I 11371 . finds largely unpolarized J /if/ 
in the helicity frame, which is in contradiction to both 
the CSM and NRQCD predictions. The early ALICE 
data E3 however favors NRQCD over the CSM. We 
look very much forward to the forthcoming, more pre- 
cise polarization measurements at the LHC, which have 
the potential to clearly confirm or dismiss the various 
suggested J/i/j production mechanisms. 

We also note that in both photo- and hadroproduc- 
tion the perturbative expansion in a s seems to converge 



much more rapidly after inclusion of the CO states as 
compared to the CSM, where sign and shape of the vari- 
ous distributions are radically changed when going from 
LOtoNLO. 

We very much appreciate that very recently a sec- 
ond, independent, NLO NRQCD analysis involving the 
polarization parameter Ag in J /if/ hadroproduction has 
been performed BTll . In that preprint it was shown that, 
neglecting data from all production mechanisms except 
hadroproduction, sets of CO LDMEs can be constructed 
that describe both the measured hadroproduction yield 
and values of Ag in the helicity frame close to zero. Un- 
fortunately, as emphasized in that preprint, even after 
including the CDF polarization data it is still not pos- 
sible to extract the three CO LDMEs independently in 
a hadroproduction-only fit, let alone to test their univer- 
sality. The analysis |41 1 does however allow us to cross- 
check parts of our results, namely dcroo and dcr u in 
the helicity frame. Fortunately, both works agree there, 
since we can reproduce their figure 2, and when using 
the LDME set in the first row of their table 1 , our pre- 
diction for Ag nicely lies within the "NLO total" band of 
their figure 1 . Furthermore, we have verified that none 
of their proposed LDME sets is compatible with the ob- 
served photoproduction yield. Thus both works agree 
on that the observed hadro- and photoproduction yields 
and the CDF polarization measurement [ 38 1 can not be 
simultaneously described by a single set of CO LDMEs. 

5. Summary 

We have presented a rigorous NLO analysis of the 
J/if/ yield and polarization in various production mech- 
anisms within the factorization theorem of NRQCD. We 
have extracted values for the three CO LDMEs by fitting 
to 194 points of inclusive J /if/ production data from dif- 
ferent hadroproduction, photoproduction, two-photon 
scattering and electron-positron annihilation measure- 
ments. We have shown that this fit is constrained and 
independent of the choice of possible lower cuts on the 
J/ijj transverse momenta. Our global fit can describe 
all data well, except perhaps the two-photon scattering 
data by DELPHI. Besides that, the extracted LDMEs do 
exhibit the scaling behavior predicted by the NRQCD 
scaling rules. 

In a second step we have used the CO LDMEs thus 
extracted to make predictions for J /if/ polarization ob- 
servables in photo- and hadroproduction and compare 
with the currently available experimental data. As for 
photoproduction, HERA data is not precise enough for 
definite conclusions, here we have to wait for future ex- 
periments, possibly at the LHeC. But as for hadropro- 
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Figure 6: Polarization parameters Ag and A^, for hadroproduction using CO LDME set B of table[T] compared to data measured by CDF at Tevatron 
in run I [ 37 ] and II [ 38 1 and by ALICE at the LHC ['391. Blue bands: Uncertainties of NLO CS curve due to scale variations. Yellow bands: 
Uncertainties of NLO CS+CO curve due to scale variations and uncertainties of the CO LDMEs. From 1 16 1. 
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Figure 7: The predictions for do-go, dcr\\ and dcr\ -\ in the helicity and Collins-Soper frames decomposed into the contributions of the individual 
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duction, the precise results for Ag from the CDF col- 
laboration measured at Tevatron run II in the helicity 
frame are in stark conflict with the NRQCD predictions. 
CDF measured unpolarized J/iJ/ here while NRQCD 
predicts a strong transverse polarization. However, the 
CDF measurements at Tevatron run I and II seem to be 
partly in conflict with each other, and the early data from 
the ALICE collaboration at low transverse momentum 
at the LHC is compatible with NRQCD. Future, more 
precise polarization measurements at the LHC, also at 
higher transverse momenta, will have the potential to 
clearly confirm or dismiss the LDME universality. 

We stress that our present analysis still ignores feed- 
down contributions from higher charmonium states. Al- 
though they have been taken care of in our global fit 
through a subtraction of estimated contributions from 
the prompt experimental data, and their influence on the 
polarization observables is expected to be very minor, a 
rigorous inclusion of the feed-down contributions is still 
left for future work. 
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